Use of synteny to identify candidate genes underlying QTL controlling stomatal traits in faba bean (Vicia faba L.)

Key message We have identified QTLs for stomatal characteristics on chromosome II of faba bean by applying SNPs derived from M. truncatula , and have identified candidate genes within these QTLs using synteny between the two species. Abstract Faba bean (Vicia faba L.) is a valuable food and feed crop worldwide, but drought often limits its production, and its genome is large and poorly mapped. No information is available on the effects of genomic regions and genes on drought adaptation characters such as stomatal characteristics in this species, but the synteny between the sequenced model legume, Medicago truncatula, and faba bean can be used to identify candidate genes. A mapping population of 211 F5 recombinant inbred lines (Mélodie/2 × ILB 938/2) were phenotyped to identify quantitative trait loci (QTL) affecting stomatal morphology and function, along with seed weight, under well-watered conditions in a climate-controlled glasshouse in 2013 and 2014. Canopy temperature (CT) was evaluated in 2013 under water-deficit (CTd). In total, 188 polymorphic single nucleotide polymorphisms (SNPs), developed from M. truncatula genome data, were assigned to nine linkage groups that covered ~928 cM of the faba bean genome with an average inter-marker distance of 5.8 cM. 15 putative QTLs were detected, of which eight (affecting stomatal density, length and conductance and CT) co-located on chromosome II, in the vicinity of a possible candidate gene—a receptor-like protein kinase found in the syntenic interval of M. truncatula chromosome IV. A ribose-phosphate pyrophosphokinase from M. truncatula chromosome V, postulated as a possible candidate gene for the QTL for CTd, was found some distance away in the same chromosome. These results demonstrate that genomic information from M. truncatula can successfully be translated to the faba bean genome.

The identification of QTL controlling ergot sclerotia size in hexaploid wheat implicates a role for the Rht dwarfing alleles

The fungal pathogen Claviceps purpurea infects ovaries of a broad range of temperate grasses and cereals, including hexaploid wheat, causing a disease commonly known as ergot. Sclerotia produced in place of seed carry a cocktail of harmful alkaloid compounds that result in a range of symptoms in humans and animals, causing ergotism. Following a field assessment of C. purpurea infection in winter wheat, two varieties ‘Robigus’ and ‘Solstice’ were selected which consistently produced the largest differential effect on ergot sclerotia weights. They were crossed to produce a doubled haploid mapping population, and a marker map, consisting of 714 genetic loci and a total length of 2895 cM was produced. Four ergot reducing QTL were identified using both sclerotia weight and size as phenotypic parameters; QCp.niab.2A and QCp.niab.4B being detected in the wheat variety ‘Robigus’, and QCp.niab.6A and QCp.niab.4D in the variety ‘Solstice’. The ergot resistance QTL QCp.niab.4B and QCp.niab.4D peaks mapped to the same markers as the known reduced height (Rht) loci on chromosomes 4B and 4D, Rht-B1 and Rht-D1, respectively. In both cases, the reduction in sclerotia weight and size was associated with the semi-dwarfing alleles, Rht-B1b from ‘Robigus’ and Rht-D1b from ‘Solstice’. Two-dimensional, two-QTL scans identified significant additive interactions between QTL QCp.niab.4B and QCp.niab.4D, and between QCp.niab.2A and QCp.niab.4B when looking at sclerotia size, but not between QCp.niab.2A and QCp.niab.4D. The two plant height QTL, QPh.niab.4B and QPh.niab.4D, which mapped to the same locations as QCp.niab.4B and QCp.niab.4D, also displayed significant genetic interactions.

Reduced height (Rht) and photoperiod insensitivity (Ppd) allele associations with establishment and early growth of wheat in contrasting production systems

Near isogenic lines (NILs) varying for genes for reduced height (Rht) and photoperiod insensitivity (Ppd-D1a) in a cv. Mercia background (rht (tall), Rht-B1b, Rht-D1b, Rht-B1c, Rht8c + Ppd-D1a, Rht-D1c, Rht12) were compared at one field site but within contrasting ('organic' vs. 'conventional') rotational and agronomic contexts, in each of 3 years. In the final year, further NILs (rht (tall), Rht-B1b, Rht-D1b, Rht-B1c, Rht-B1b + Rht-D1b, Rht-D1b + Rht-B1c) in both Maris Huntsman and Maris Widgeon backgrounds were added together with 64 lines of a doubled haploid (DH) population [Savannah (Rht-D1b) x Renesansa (Rht-8c + Ppd-D1a)]. Assessments included laboratory tests of germination and coleoptile length, and various field measurements of crop growth between emergence and pre jointing [plant population, tillering, leaf length, ground cover (GC), interception of photosynthetically active radiation (PAR), crop dry matter (DM) and nitrogen accumulation (N), far red: red reflectance ratio (FR:R), crop height, and weed dry matter]. All of the dwarfing alleles except Rht12 in the Mercia background and Rht8c in the DHs were associated with reduced coleoptile length. Most of the dwarfing alleles (depending on background) reduced seed viability. Severe dwarfing alleles (Rht-B1c, Rht-D1c and Rht12) were routinely associated with fewer plant numbers and reduced early crop growth (GC, PAR, DM, N, FR:R), and in 1 year, increased weed DM. In the Mercia background and the DHs the semi-dwarfing allele Rht-D1b was also sometimes associated with reductions in early crop growth; no such negative effects were associated with the marker for Rht8c. When significant interactions between cropping system and genotype did occur it was because differences between lines were more exaggerated in the organic system than in the conventional system. Ppd-D1a was associated positively with plant numbers surviving the winter and early crop growth (GC, FR:R, DM, N, PAR, height), and was the most significant locus in a QTL analysis. We conclude that, within these environmental and system contexts, genes moderating development are likely to be more important in influencing early resource capture than using Rht8c as an alternative semi-dwarfing gene to Rht-D1b.

A novel variant on chromosome 7q22.3 associated with mean platelet volume, counts, and function

Mean platelet volume (MPV) and platelet count (PLT) are highly heritable and tightly regulated traits. We performed a genome-wide association study for MPV and identified one SNP, rs342293, as having highly significant and reproducible association with MPV (per-G allele effect 0.016 +/- 0.001 log fL; P < 1.08 x 10(-24)) and PLT (per-G effect -4.55 +/- 0.80 10(9)/L; P < 7.19 x 10(-8)) in 8586 healthy subjects. Whole-genome expression analysis in the 1-MB region showed a significant association with platelet transcript levels for PIK3CG (n = 35; P = .047). The G allele at rs342293 was also associated with decreased binding of annexin V to platelets activated with collagen-related peptide (n = 84; P = .003). The region 7q22.3 identifies the first QTL influencing platelet volume, counts, and function in healthy subjects. Notably, the association signal maps to a chromosome region implicated in myeloid malignancies, indicating this site as an important regulatory site for hematopoiesis. The identification of loci regulating MPV by this and other studies will increase our insight in the processes of megakaryopoiesis and proplatelet formation, and it may aid the identification of genes that are somatically mutated in essential thrombocytosis. (Blood. 2009; 113: 3831-3837)

Mapping quantitative trait loci for flag leave senescence as a yield determinant in winter wheat under optimal and drought stressed environments

The timing of flag leaf senescence (FLS) is an important determinant of yield under stress and optimal environments. A doubled haploid population derived from crossing the photo period-sensitive variety Beaver,with the photo period-insensitive variety Soissons, varied significantly for this trait, measured as the percent green flag leaf area remaining at 14 days and 35 days after anthesis. This trait also showed a significantly positive correlation with yield under variable environmental regimes. QTL analysis based on a genetic map derived from 48 doubled haploid lines using amplified fragment length polymorphism (AFLP) and simple sequence repeat (SSR) markers, revealed the genetic control of this trait. The coincidence of QTL for senescence on chromosomes 2B and 2D under drought-stressed and optimal environments, respectively, indicate a complex genetic mechanism of this trait involving the re-mobilisation of resources from the source to the sink during senescence.

Genetic analysis of heat tolerance at anthesis in rice

Genetic analysis of heat tolerance will help breeders produce rice (Oryza sativa L.) varieties adapted to future climates. An F6 population of 181 recombinant inbred lines of Bala (tolerant) × Azucena (susceptible) was screened for heat tolerance at anthesis by measuring spikelet fertility at 30°C (control) and 38°C (high temperature) in experiments conducted in the Philippines and the United Kingdom. The parents varied significantly for absolute spikelet fertility under control (79–87%) and at high temperature (2.9–47.1%), and for relative spikelet fertility (high temperature/control) at high temperature (3.7–54.9%). There was no correlation between spikelet fertility in control and high-temperature conditions and no common quantitative trait loci (QTLs) were identified. Two QTLs for spikelet fertility under control conditions were identified on chromosomes 2 and 4. Eight QTLs for spikelet fertility under high-temperature conditions were identified on chromosomes 1, 2, 3, 8, 10, and 11. The most significant heat-responsive QTL, contributed by Bala and explaining up to 18% of the phenotypic variation, was identified on chromosome 1 (38.35 mega base pairs on the rice physical genome map). This QTL was also found to influence plant height, explaining 36.6% of the phenotypic variation. A comparison with other studies of abiotic (drought, cold, salinity) stresses showed QTLs at similar positions on chromosomes 1, 3, 8, and 10, suggesting common underlying stress-responsive regions of the genome.

Adaptation of tree growth to elevated CO2: quantitative trait loci for biomass in Populus

Information on the genetic variation of plant response to elevated CO2 (e[CO2]) is needed to understand plant adaptation and to pinpoint likely evolutionary response to future high atmospheric CO2 concentrations.• Here, quantitative trait loci (QTL) for above- and below-ground tree growth were determined in a pedigree – an F2 hybrid of poplar (Populus trichocarpa and Populus deltoides), following season-long exposure to either current day ambient CO2 (a[CO2]) or e[CO2] at 600 µl l−1, and genotype by environment interactions investigated.• In the F2 generation, both above- and below-ground growth showed a significant increase in e[CO2]. Three areas of the genome on linkage groups I, IX and XII were identified as important in determining above-ground growth response to e[CO2], while an additional three areas of the genome on linkage groups IV, XVI and XIX appeared important in determining root growth response to e[CO2].• These results quantify and identify genetic variation in response to e[CO2] and provide an insight into genomic response to the changing environment

Development of a dense SNP-based linkage map of an apple rootstock progeny using the Malus Infinium whole genome genotyping array

Background A whole-genome genotyping array has previously been developed for Malus using SNP data from 28 Malus genotypes. This array offers the prospect of high throughput genotyping and linkage map development for any given Malus progeny. To test the applicability of the array for mapping in diverse Malus genotypes, we applied the array to the construction of a SNPbased linkage map of an apple rootstock progeny. Results Of the 7,867 Malus SNP markers on the array, 1,823 (23.2 %) were heterozygous in one of the two parents of the progeny, 1,007 (12.8 %) were heterozygous in both parental genotypes, whilst just 2.8 % of the 921 Pyrus SNPs were heterozygous. A linkage map spanning 1,282.2 cM was produced comprising 2,272 SNP markers, 306 SSR markers and the S-locus. The length of the M432 linkage map was increased by 52.7 cM with the addition of the SNP markers, whilst marker density increased from 3.8 cM/marker to 0.5 cM/marker. Just three regions in excess of 10 cM remain where no markers were mapped. We compared the positions of the mapped SNP markers on the M432 map with their predicted positions on the ‘Golden Delicious’ genome sequence. A total of 311 markers (13.7 % of all mapped markers) mapped to positions that conflicted with their predicted positions on the ‘Golden Delicious’ pseudo-chromosomes, indicating the presence of paralogous genomic regions or misassignments of genome sequence contigs during the assembly and anchoring of the genome sequence. Conclusions We incorporated data for the 2,272 SNP markers onto the map of the M432 progeny and have presented the most complete and saturated map of the full 17 linkage groups of M. pumila to date. The data were generated rapidly in a high-throughput semi-automated pipeline, permitting significant savings in time and cost over linkage map construction using microsatellites. The application of the array will permit linkage maps to be developed for QTL analyses in a cost-effective manner, and the identification of SNPs that have been assigned erroneous positions on the ‘Golden Delicious’ reference sequence will assist in the continued improvement of the genome sequence assembly for that variety.

Segmental chromosomal duplications harbouring group IV CONSTANS-like genes in cereals

Comparative mapping is an important component of map-based cloning in large-genome cereal species. We describe evidence of a segmental chromosomal duplication harbouring CONSTANS-like genes in barley that predates the divergence of the Oryzoideae (rice) and Pooideae (brachypodium, barley, wheat) clades, and discuss the implications of such events for comparative mapping and QTL cloning in temperate cereal crops.

Molecular barley breeding

Breeding progress in barley yield in the UK is being sustained at a rate in the order of 1% per annum against a background of declining seed sales. Commercial barley breeders are largely concentrating upon the elite local gene pool but with genotypic evidence suggesting that there is still considerable variation between current recommended cultivars, even those produced as half-sibs by the same breeder. Marker Assisted Selection (MAS) protocols could be substituted for conventional selection for a number of major-gene targets but, in the majority of cases, conventional selection is more resource efficient. Results from current QTL mapping studies have not yet identified sufficiently robust and validated targets for UK barley breeders to adopt MAS to assist in the selection of complex traits such as yield and malting quality. Results from multiple population mapping amongst the elite gene pool being utilised by breeders and from association studies of elite germplasm tested as part of the UK recommended list trial process do, however, show some promise.

Shoot calcium and magnesium concentrations differ between subtaxa, are highly heritable, and associate with potentially pleiotropic loci in Brassica oleracea

Calcium (Ca) and magnesium (Mg) are the most abundant group II elements in both plants and animals. Genetic variation in shoot Ca and shoot Mg concentration (shoot Ca and Mg) in plants can be exploited to biofortify food crops and thereby increase dietary Ca and Mg intake for humans and livestock. We present a comprehensive analysis of within-species genetic variation for shoot Ca and Mg, demonstrating that shoot mineral concentration differs significantly between subtaxa (varietas). We established a structured diversity foundation set of 376 accessions to capture a high proportion of species-wide allelic diversity within domesticated Brassica oleracea, including representation of wild relatives (C genome, 1n = 9) from natural populations. These accessions and 74 modern F-1 hybrid cultivars were grown in glasshouse and field environments. Shoot Ca and Mg varied 2- and 2.3-fold, respectively, and was typically not inversely correlated with shoot biomass, within most subtaxa. The closely related capitata (cabbage) and sabauda (Savoy cabbage) subtaxa consistently had the highest mean shoot Ca and Mg. Shoot Ca and Mg in glasshouse-grown plants was highly correlated with data from the field. To understand and dissect the genetic basis of variation in shoot Ca and Mg, we studied homozygous lines from a segregating B. oleracea mapping population. Shoot Ca and Mg was highly heritable (up to 40). Quantitative trait loci (QTL) for shoot Ca and Mg were detected on chromosomes C2, C6, C7, C8, and, in particular, C9, where QTL accounted for 14 to 55 of the total genetic variance. The presence of QTL on C9 was substantiated by scoring recurrent backcross substitution lines, derived from the same parents. This also greatly increased the map resolution, with strong evidence that a 4-cM region on C9 influences shoot Ca. This region corresponds to a 0.41-Mb region on Arabidopsis (Arabidopsis thaliana) chromosome 5 that includes 106 genes. There is also evidence that pleiotropic loci on C8 and C9 affect shoot Ca and Mg. Map-based cloning of these loci will reveal how shoot-level phenotypes relate to Ca 21 and Mg 21 uptake and homeostasis at the molecular level.

Shoot zinc (Zn) concentration varies widely within Brassica oleracea L. and is affected by soil Zn and phosphorus (P) levels

The low availability of zinc (Zn) in soils and crops affects dietary Zn intake worldwide. This study sought to determine if the natural genetic variation in shoot Zn concentrations (Zn(shoot)) is sufficient to pursue a crop improvement breeding strategy in a leafy vegetable crop. The gene-pool of Brassica oleracea L. was sampled using a large (n = 376) diversity foundation set (DFS), representing almost all species-wide common allelic variation, and 74 commercial varieties (mostly F(1)). The DFS genotypes were grown at low and high soil phosphorus (P) levels under glasshouse and field conditions, and also in a Zn-deficient soil, with or without Zn-fertilisation, in a glasshouse. Despite the large variation in Zn(shoot) among genotypes, environment had a profound effect on Zn(shoot) The heritability of Zn(shoot) was significant, but relatively low, among 90 doubled-haploid (DH) lines from a mapping population. While several quantitative trait loci (QTL) associated with Zn(shoot) occurred on chromosomes C2, C3, C5, C7, and C9, these were generally weak and conditional upon growth conditions. Breeding for Zn(shoot) in B. oleracea is therefore likely to be challenging. Shoot P concentrations increased substantially in all genotypes under low soil Zn conditions. Conversely, only some genotypes had increased Zn(shoot) at low soil P levels. Sufficient natural genetic variation may therefore exist to study some of the interactions between Zn and P nutrition.

Regulatory hotspots are associated with plant gene expression under varying soil phosphorus supply in Brassica rapa

Gene expression is a quantitative trait that can be mapped genetically in structured populations to identify expression quantitative trait loci (eQTL). Genes and regulatory networks underlying complex traits can subsequently be inferred. Using a recently released genome sequence, we have defined cis- and trans-eQTL and their environmental response to low phosphorus (P) availability within a complex plant genome and found hotspots of trans-eQTL within the genome. Interval mapping, using P supply as a covariate, revealed 18,876 eQTL. trans-eQTL hotspots occurred on chromosomes A06 and A01 within Brassica rapa; these were enriched with P metabolism-related Gene Ontology terms (A06) as well as chloroplast-and photosynthesis-related terms (A01). We have also attributed heritability components to measures of gene expression across environments, allowing the identification of novel gene expression markers and gene expression changes associated with low P availability. Informative gene expression markers were used to map eQTL and P use efficiency-related QTL. Genes responsive to P supply had large environmental and heritable variance components. Regulatory loci and genes associated with P use efficiency identified through eQTL analysis are potential targets for further characterization and may have potential for crop improvement.

Natural genetic variation in caesium (Cs) accumulation by Arabidopsis thaliana

Ingestion of caesium (Cs) radioisotopes poses a health risk to humans. Crop varieties that accumulate less Cs in their edible tissues may provide a useful countermeasure. This study was performed to determine whether quantitative genetics on a model plant (Arabidopsis thaliana) might inform such 'safe'-crop strategies. Arabidopsis accessions and recombinant inbred lines (RILs), from Landsberg erecta (Ler) x Cape Verdi Island (Cvi), Ler x Columbia (Col), and Niederzenz (Nd) x Col mapping populations, were grown on agar supplemented with subtoxic levels of Cs. Shoot Cs concentration varied up to three-fold, and shoot f. wt varied up to 25-fold within populations. The heritability of growth and Cs accumulation traits ranged from 0.06 to 0.28. Four quantitative trait loci (QTL) accounted for > 80 of the genetic contribution to the total phenotypic variation in shoot Cs concentration in the Ler x Col population. QTL identified in this study, in particular, QTL co-localizing to the top and bottom regions of Chromosomes I and V in two different mapping populations, are amenable to positional cloning and, through collinearity, may inform selection or breeding strategies for the development of 'safe' crops.

High-throughput root phenotyping screens identify genetic loci associated with root architectural traits in Brassica napus under contrasting phosphate availabilities

Background and Aims: Phosphate (Pi) deficiency in soils is a major limiting factor for crop growth worldwide. Plant growth under low Pi conditions correlates with root architectural traits and it may therefore be possible to select these traits for crop improvement. The aim of this study was to characterize root architectural traits, and to test quantitative trait loci (QTL) associated with these traits, under low Pi (LP) and high Pi (HP) availability in Brassica napus. Methods: Root architectural traits were characterized in seedlings of a double haploid (DH) mapping population (n = 190) of B. napus 'Tapidor' x 'Ningyou 7' (TNDH) using high-throughput phenotyping methods. Primary root length (PRL), lateral root length (LRL), lateral root number (LRN), lateral root density (LRD) and biomass traits were measured 12 d post-germination in agar at LP and HP. Key Results: In general, root and biomass traits were highly correlated under LP and HP conditions. 'Ningyou 7' had greater LRL, LRN and LRD than 'Tapidor', at both LP and HP availability, but smaller PRL. A cluster of highly significant QTL for LRN, LRD and biomass traits at LP availability were identified on chromosome A03; QTL for PRL were identified on chromosomes A07 and C06. Conclusions: High-throughput phenotyping of Brassica can be used to identify root architectural traits which correlate with shoot biomass. It is feasible that these traits could be used in crop improvement strategies. The identification of QTL linked to root traits under LP and HP conditions provides further insights on the genetic basis of plant tolerance to P deficiency, and these QTL warrant further dissection.